A1 = readmatrix('NiI2_20nm_20V.csv')
A2 = readmatrix('NiI2_20nm_-20V.csv')

x1 = A1(:,1)
y1 = A1(:,2)*1e15

x2 = A2(:,1)
y2 = A2(:,2)*1e15

hold on
xlim([50, 65])
plot(x1,y1,'Linewidth',1,'Color',[0.8 0 0])
plot(x2,y2,'Linewidth',1,'Color',[0 0 0.8])
xlabel('Temperature [K]')
ylabel('Pyroelectric Current [fA]')
ylim([10,70])
legend('20 kV/m', '-20 kV/m')
set(gca,'XMinorTick','on','YMinorTick','on')
legend boxoff
box on
hold off

set(gcf,'Position',[0, 0, 200, 200])